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Abstract 

We consider the proposal that dark matter (DM) is composed of a spin-3/2 particle that is a singlet of 
the standard model (SM). Its leading effective interactions with ordinary matter involve a pair of their 
fields and a pair of SM fermions, in the form of products of chiral currents. We make a comprehensive 
analysis on possible phenomenological effects of the interactions in various experiments and observations. 
These include collider searches for monojet plus missing transverse energy events, direct detections of DM 
scattering off nuclei, possible impacts on the gamma rays and antiproton-to-proton flux ratio in cosmic rays, 
and the obsei^ved relic density. The current data akeady set strong constraints on the effective interactions 
in a complementary manner. The constraint from collider searches is most effective at a relatively low 
mass of DM, and the antiproton-to-proton flux ratio offers the best bound for a heavy DM, while the spin- 
independent direct detection is the best in between. For DM mass of order 10 GeV to 1 TeV, the effective 
interaction scale is constrained to be typically above a few tens TeV. 



*Electronic address: dingran@mail.nankai.edu.cn 



^Electronic address: liaoy@nankai.edu.cn 



^Electronic address: liujy@tjut.edu.cn 



^Electronic address: |wangkail @zju.edu"xn 



I. INTRODUCTION 



The evidence for the domination of dark matter (DM) over ordinary matter in our universe is 
robust, but still restricted to its gravitational effects after years of efforts, from Zwicky's sugges- 
tion in 1930's to explain the rotation curves of galaxies and galaxy clusters to recent precision 
measurements on cosmic microwave background; see [111 for brief reviews. If some part of DM 
is nonbaryonic as data indicated and has the nature of particles, it should interact weakly with 
ordinary matter to cause other effects that could be observable by the means of particle detection. 
Indeed, there have been many observational and experimental activities trying to reveal various 
aspects of DM particles, from direct and indirect detections to collider searches. They have pro- 
vided useful constraints on the nature of DM particles, and may hopefully discover them in the 
near future. 

Many new physics models contain massive neutral particles whose stability is protected by cer- 
tain exact or approximate symmetries, and thus could serve as DM particles. Most extensively 
studied are perhaps supersymmetric models; also popular are models based on extra dimensions 
yD-da], and little-Higgs models [|7|, |8D, to mention a few among many; see Il9l1- [|12n for detailed 
reviews. Since the basic properties of DM particles are more or less fixed in these models, for 
instance, their spins, structures and orders of magnitude of interactions with ordinary matter, it is 
possible to make rather detailed predictions on their observational effects. On the other hand, the 
physical relevance of the models themselves remains to be experimentally verified. Considering 
our still limited knowledge on DM particles, it is necessary to avoid theoretical biases in exploring 
various possibilities. In such a circumstance, the effective field theory approach could be very 
useful lll3n . By assuming basic properties of a DM particle such as its spin and mass, one exhausts 
its effective interactions with ordinary matter that respect known symmetries and are of leading 
order at low energies while leaving interaction strengths as phenomenological parameters. The 
physical effects can then be determined in terms of those parameters and confronted with experi- 
mental measurements. If a DM particle is fortunately discovered, the rough information gathered 
for those properties and parameters could be employed as important physical input in planning 
future facilities to reveal its underlying dynamics. 

The DM candidates of a spin zero, spin- 1/2, and spin one particle have been exhaustively 
studied in the literature in the framework of effective field theory lll4ll - [|30n . Recently, two of 
us have considered the possibility that the DM particle may have spin-3/2 Bill (see also Ref. 
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13211 on the same suggestion). In that work, the quantum numbers of the DM particle under the 
standard model (SM) gauge group were not specified, and the constraints from direct and indirect 
detection data were found similar to those for a spin- 1/2 particle. Since DM interacts very weakly 
with ordinary matter, it should more naturally be a SM singlet. Here we consider this option 
and investigate the constraints coming from collider measurements as well as direct and indirect 



detections. A light, singlet DM particle of spin-3/2 has also been studied earlier [ |33( | in rare 
decays of the K and B mesons where the particle appears as missing energy in final states. While 
the kinematics of such a particle is similar to that of a gravitino, which is also a DM candidate in 
supergravity models (see [34] for a review), the interactions to be examined here are very different. 



A charged spin-3/2 particle was also proposed earlier [13511 as a constituent of the so-called dark 



atoms. More recently, a specific model of a spin-3/2 particle was suggested and its direct detection 



examined 113611 . in which the particle is charged under the SM gauge group. The possible relevance 
of spin-3/2 particles has also been considered in other contexts, see for instance, Ref. [jj/l, on 
collider effects of a spin-3/2 top partner. 

The paper is organized as follows. In the next section we considerpossible effective interactions 
of a singlet, spin-3/2 particle with the standard model particles, and set up our conventions for 
spin-3/2 particles. This is followed by sec Union the Large Hadron Collider (LHC) effects of DM 
particles that may appear as missing energy in monojet events. In sections |IV] and |Vl we consider, 
respectively, the direct detection via DM scattering off nuclei and the indirect detection through 
impacts on the cosmic rays. All of these constraints are combined in sec |VI] together with that 
from the observed relic density. We summarize briefly in the last section. 

II. EFFECTIVE INTERACTIONS 

The field corresponding to a particle of spin-3/2 and mass M is described by a vector- 



spinor, with the constraint, y^^^ = 13811 . The free field satisfies the equation of mo 



tion, (i^ — M)^^ = 0. The wavefunction of such a particle with momentum p and helicity A, 
U^{p,X), can be constructed from that of a Dirac particle and the polarization of a spin one 



particle in terms of the Clebsch-Gordan coefficients [|39l1 . For later applications we will need 
the polarization sums for such a particle, Pv^{p) = Ea Uv{p-,^)U^{p,X), and for an antiparticle, 
Qviiip) = 'Lxyv{P:X)y^{PiX), with V^{p,X) being the antiparticle's wavefunction. They are 
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known to be 



Piivip) = -(^+M)(^T^y{p)-^fTp^{p)TyMfy (1) 

and Q^vip) = P^v{p)\m^-m, where T^v{p) = giiv - PiiPv / P^ and p^=M^. 

We consider the effective interactions of a spin-3/2 particle with the SM fermions, 

U (-1), Er (-2); Qi^ (1/3), Ur (4/3), Dr (-2/3); (2) 

where the number in parentheses denotes the hypercharge Y which is related to the electric charge 
Q and third weak i so spin hy the convention Q = + Y /2. We start with the operators that 
involve a pair of SM fermions and a pair of spin-3/2 fields. Lorentz invariance allows for a list 
of fourteen independent structures [13 ill as explicitly verified using the generalized Fierz identities 



4011 ■ Demanding to be a SM singlet reduces the list to the following four: 



'!'^fP+'i"'fLrcfL, (3) 



where P± = {1 iy^) /2 and (/r) refers to any SM doublet (singlet) fermion field. Note that the 
operators are automatically flavor diagonal and that reshuffling a with an / does not introduce 



independent operators according to [HOll . The corresponding effective interactions are parameter- 
ized as 

^eff = (c« + c{^{)+A-2 £ (c{4 + c{4), (4) 

f=L,Q f=E,U,D 

where A is the typical energy scale inducing the interactions and q s are dimensionless real pa- 
rameters presumably of order one. Our later numerical analysis will be based on this effective 
Lagrangian. To reduce the number of unknowns, we follow the usual practice: we treat one oper- 
ator at a time and assume a universal for all relevant SM fermions. 

We mention briefly some other operators that can be built out of the spin-3/2 field and the SM 
fields using the approach in [|40|l . Each of these operators violates either the lepton or baryon num- 
ber but not both, and is also forbidden if the DM particle carries certain conserved parity. They 
could thus be phenomenologically dangerous, and we will study them elsewhere. An operator 



involving a single SM fermion requires an odd number of the field. With a single such 
an operator has the lowest possible dimension five, {D^H)^^'^Li, plus one constructed with the 
help of charge conjugation. Here H is the Higgs doublet with H = eH*, and is the SM gauge 
covariant derivative. There is no similar operator involving a quark field. It is not possible either 
to form a dimension six operator involving a single SM fermion and three ^^s without including 
a genuinely neutral field of neutrinos, V^. If Vr is indeed introduced, there are then seven such op- 



erators according to the results in [|40ll . An operator involving three SM fermions contains at least 
a single corresponding to a dimension six operator. The situation is a bit complicated since 
one can have pure-lepton, pure-quark, and mixed lepton-quark operators, involving additional 
color contraction for the latter two. For simplicity, we show here only the pure-lepton operators. 
Lorentz invariance allows a complete and independent list of four chirality-diagonal (like 2^ 
structures and four chirality-flipped (like ^{^) ones. But gauge symmetry of SM singles out only 



the following operators, e"''Lfj^o^''EkRL^j^rv^^, e^'^iLf^jCo^^L^j^EkRYv^^, where /, j, k stand 
for family and a, b for the third weak isospin. 

III. DIRECT DM PRODUCTION AT LHC AND ITS CONSTRAINTS 

In this section, we study the collider phenomenology of spin-3/2 DM using the effective in- 
teractions shown in sec III Since the DM particle is electrically and chromatically neutral, it is 
completely invisible for detectors and only appears as missing transverse energy (Xr) at LHC. 
The direct production of DM pairs would then be completely invisible with nothing to trigger on. 
For the trigger purpose, we could focus on the production of a DM pair in association with an ini- 
tial state radiation jet or photon. In the case of a monojet plus final state, the level one trigger 
requires that the sum of the jet transverse momentum {pj) and missing transverse energy, pj +Xr, 
be greater than 250 GeV or so at the LHC detectors. Studies on monojet signatures with spin- 1/2 
DM effective operators have been performed extensively Jm, 41-47]. Here, we use the latest data 
on search of monojet final states to constrain the effective interactions involving spin-3/2 DM. 

There are three independent subprocesses for production of a monojet plus a DM pair, 

qq^g^'^'iy and gq{q) ^ q{q)^' ^ , (5) 

corresponding to Feynman diagrams in Fig. [T] where g stands for a gluon. We compute in the 
appendix the spin- and color-summed and -averaged matrix elements squared due to various op- 
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FIG. 1 : Feynman diagrams for a monojet plus a spin-3/2 DM pair production at a hadron collider. 



erators To illustrate the feature, we show here the result for the subprocess q{j>\)q{j>2) — >■ 
gikjY^ ^{k\fVy{k2) due to the operator G\, i.e., assuming c\ = \ and C2,3,4 = 0, 

„2 o 

y 1^12 = ii_ — ^ — B. (6) 

Here gs is the QCD gauge coupling, and assuming massless quarks, one has, 

B = [4M^(l-xi)2 + 2M45(l-;ci-X2)(l-JCi-yi)(5-5jci-7yi) 

-4MV(1 -XX -X2f[\ -XI -yif + s\\ -xi -X2)\\ -xi -yif 
+{^\^X2,yi^y2), (7) 

where s = {p\+ piY, and x,, are kinematical variables defined in the appendix. The result due 
to the operator ^2 alone is identical. The same degeneracy also occurs between the operators 
and ^4, and is due to the spin summation and averaging. Our numerical analysis shows that the 
difference in angular distributions between the operators ^1 2 on one side and 1^3 4 on the other, as 
the expressions for various 5s indicate, is further smeared out by phase space integration, resulting 
in the same total cross sections within about one percent. The four operators are thus equivalent 
for the phenomenology considered here when a universal coupling c{ is assumed, and from now 
on we will denote them simply by Compared to the case of spin- 1/2 DM, the production of 
spin-3/2 DM is enhanced in the low mass or high energy region, and we thus expect more stringent 
bounds could be set in this case. 

In Fig. [2] we show the numerical result for the total cross section of a monojet plus Xr due to 
any one operator ^ at 7 TeV LHC using different jet pj cuts. The effective scale A is fixed at 



6 



10'' 



10" 



10"' 



10"' 



X3 

b 10-9 



10 



10 



12 



15 



10" 



10 



LHC, Vs = 7TeV 
A= lOTeV 



50 100 



Pj > 220 GeV, I ri I < 2 




500 1000 3000 



M (GeV) 



FIG. 2: Total partonic cross section of direct production of a DM pair plus a monojet at 7 TeV LHC as a 
function of DM mass for any one operator iff. We choose the PDF set CTEQ6L1 with renomialization scale 
Hr = Q and factorization scale iXp = 2/2, and A = 10 TeV. 

n 

10 TeV, and we choose the parton distribution function (PDF) set CTEQ6L1 [48], with renormal- 
ization scale Hr = Q and factorization scale jlf = Q/2. The monojet cut is within |77|y < 2 with 
p^j- > 20 GeV and pj > 220 GeV respectively. The cross section is dominated in the high energy 
region by the terms in the amplitude squared; only when the DM mass becomes non-negligible 
compared with ^/s, does the phase space suppression begin to play a significant role and reduce 
the cross section rapidly. This is reminiscent of the well-known gravitino-goldstino equivalence 
in supergravity models [49I-I5II. In the high energy limit, the polarization sum Pjxv{p) approaches 
~^gjiv +2/(3M^);^p^pv, where the first and second term can be identified with the helicity states 
A = ±3/2 and A = ±1/2 respectively [|52|] . The dominance of the latter amounts to an effective 



54|] . The results using the full 



description by a spin-1/2 field ^ via ^'^ ^flj3M-^d^^' [|53, 
field and its effective description are compared in Fig. |3] where both the cross sections and 
their ratio are shown. It is clear that the ratio approaches unity in the low M region and drops when 
M increases beyond about 500 GeV. 

n n 

Recently, both ATLAS [|55|] and CMS |56l ] collaborations have released their studies on monojet 



plus Xr events based on data at -y/i = 7 TeV and with an integrated luminosity of 4.7 fb ^ and 
5.0 fb^^, respectively. We use these latest data to derive bounds on our effective operators. Since 
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FIG. 3: Total partonic cross section for production of a DM pair plus monojet at 7 TeV LHC is shown as 
a function of DM mass for operator using the full field or its effective description, a/2/3M^^5^*P. 
Also shown is the ratio of results in the two approaches. Same physical input as in Fig. |2] 

the monojet pj cut is typically required to be harder than 100 GeV, we simulate our signal events 
with a partonic monojet plus a DM pair without parton showers. In addition, in order to simulate 
the detector performance at the ATLAS and CMS detectors, we smear jets preformed according to 



the energy resolution [|57 



AE, 



0.8 



Ej ^Ej/GtN 
AEj 1.0 



)0.15 for ATLAS, 
)0.05 for CMS . 



(8) 



Ej y^Ej/GcY 

The SM background is taken from the ATLAS/CMS analysis including the corresponding uncer- 
tainties. The QCD jet production in principle can contribute to monojet plus final states due 
to the jet energy resolution. However, in this case, the distribution da/dM^j drops rapidly before 
Xr < 100 GeV. The leading SM background then consists of a monojet plus a Z boson with invis- 
ible Z decays or a monojet plus a boson decaying into soft leptons. Therefore, a large with 
a high pt monojet typically works as a good selection cut, and as we mentioned, the level one 
trigger is, p{ +Zt > 250 GeV. The AT 



cuts, namely, SR1/SR2/SR3/SR4 (SSl, 
data from ATLAS and CMS. 



AS/CMS searches have assumed four sets of selection 



56h . In Tabled we summarize the selection cuts and latest 
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ATLAS7TeV,4.7fb"^ 


CMS7TeV, S.Ofb"^ 


Signal region 


SRI /SR2 /SR3 /SR4 
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0T (GeV) > 


120 /220 /350 /500 


250 /300 /350 /400 


(GeV) > 


120 /220 /350 /500 


110 
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2.4 


NsM 


124000 /8800 /750 /83 


7842 /2757 /1225 /573 


<7SM 


4000 7400 /60 /14 


367 /167 /lOl /65 


A^obs 


124703 /8631 /785 /77 


7584 /2774 /1 142 /522 



TABLE I: Crucial cuts and data in the ATLAS 15511 and CMS tSw monojet plus Xt analyses. 



To obtain the collider bounds, we follow the definition in Ref. [|47l1 

2^ (A^(A,M)+A^SM-A^obsf 



(9) 



where Nsm is the number of SM background events with the uncertainty Osm covering both sta- 
tistical and systematical uncertainties, A'obs is the number of observed events, and N{A, M) is the 
event number of DM contribution from an effective operator. Here we require = 2-71 to derive 
the 90% CL exclusion bounds. Fig. |4] shows the results on the effective scale A corresponding 
to the signal regions from SRI to SR4, respectively. We note that the plots have a similar shape, 
namely, exhibit a nearly M^^ slope in the light DM mass region and begin to fall off around 
M > 500 GeV. 



IV, DIRECT DETECTION CONSTRAINTS 



The direct detection of DM measures the event rate and energy deposit in the collision of target 
nuclei (A'^) by DM particles (^^) in the local halo. It is customary to present the results in terms 
of the spin-summed and -averaged cross section for the collision at zero momentum transfer: 

1 ^ 



\6n{M + mf^] 



(10) 



spins 



where £/ is the scattering amplitude. Nevertheless, the procedure from the 'microscopic' inter- 
actions in Eq. ^ to the 'macroscopic' scattering amplitude is nontrivial. For a nice review on 
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FIG. 4: Lower bounds on A set by monojet plus Xr data from ATLAS ||55|] (left panel) and CMS (right) 



experiments. The cui^ves correspond to the signal regions from SRI to SR4 in tablelH 

the issue with potential uncertainties incurred, see Ref. [|90; for a summary of the procedure and 
a detailed analysis of the operators relevant to our discussion here, see Ref. [31]. Since the DM 
particles are nonrelativistic, they only feel the mass and spin of a nucleus. The collision can thus 
be classified into the spin-independent (SI) and spin-dependent (SD) ones. Our operators in Eq. 
^ are linear compositions of the 9th to 12th operators in [13 ill . In the nonrelativistic limit, these 
operators are dominated by the mass-mass and spin-spin terms, i.e., by the structures of 'S>Ya 
and y*^ 75 ® 7a 75, with the parity-mixed terms safely ignorable. They are thus equivalent, and each 
contributes simultaneously to the SI and SD cross sections: 

'^{bN)\ (11) 

'2 o20 



^0 



^SD 



l~L 2 
K 5 



(12) 



with jU = misM/ (m^/ +M) being the reduced mass of the ^-A'^ system. Here the effective coupling 
essentially accounts for the contributions of valence quarks of nucleons in a nucleus of mass 
number A and charge Z: 



bN = Zbp + {A-Z)b„, Z7p = A-2(2c" + c^), bn = A-'{c" + 2c''). 



—2 / „u 



(13) 



The effective coupling gi\i can also be decomposed into a sum of contributions in a nucleus of total 
spin Ji^: 



(14) 
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where A^*-"-* is the fraction of the proton (neutron) spin carried by the quark q [|59ll . (5p) and {Sn) 
are the expectation values of the total spin of protons and neutrons. When we consider the DM- 
proton (neutron) cross section, we have (Sp^,,)) = 1/2 and Jn = 1/2. For /S^^"\ we use the values 
in Ref. L180: A^? = = 0.78 ± 0.02, A^ = A'j = -0.48 ± 0.02, and Af = A^' = -0. 15 ± 0.02. 
Combining Eqs. (II 11121) . we get for any one of the operators in Eq. (|3]), 

0i = ^[{bNf + ^^JN{JN+l){gNf]. (15) 



There exist tensions among current experimental results. Both DAM A [|60ll and CoGENT 116 If] 
experiments claimed to have observed a positive signal, consistent with a DM particle of mass 
10 GeV and spin-independent cross section asi ~ 2 x 10^^^ cm^ and asi ~ 7 x 10^"^^ cm^, respec- 



tively. On the contrary, XENON [|62 



63|], CDMS m 



6511 and other experiments reported negative 



results. We do not consider the DAMA and CoGENT results in this paper, because one would need 
some isospin-violating DM model to compromise with other experiments in a consistent way. 
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FIG. 5: 90% CL upper limits inferred by ATLAS and CMS limits on SI DM-nucleon cross section as 
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a functionjof m_ass M. Also shown are 90% CL limits from the XENONIOO ^M, XENONIO IfiJ, and 
CDMSII 



6511 experiments. 



We present upper limits inferred by the ATLAS and CMS 90% CL limits on the SI and SD cross 
sections in Figs. [5] and [6] respectively. For comparison, the exclusion curves from various direct 
detection experiments are also depicted; see the captions for the detail. As one can see from the 
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figures, SD detections generally set a weaker bound than the LHC searches, while SI detections 
can yield stronger constraints than the LHC searches for a relatively large DM mass, similarly to 
the case of a spin- 1/2 DM particle. However, for a light DM particle of spin-3/2, the LHC already 
sets a much more stringent bound than the direct detections. This is in sharp contrast to the case 
of a spin- 1/2 DM particle, and originates from the enhancement discussed in sec Hill 
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FIG. 6: 90% upper limits inferred by ATLAS and CMS limits on SD DM-nucleon scattering as a func 



tion of mass M. Also shown are 90% CL limits from the SIMPLE ||66|], COUPP |l67|l, Picasso |l68|], and 



XENON 10 1691] experiments. 



V. INDIRECT DETECTION CONSTRAINTS 



When the DM particles in our Galaxy annihilate, they produce leptons and quarks which in- 
teract further with the interstellar medium to initiate more secondary particles, including gamma 
rays, neutrinos, positrons, antiprotons, etc. This provides an additional source of cosmic rays on 
top of the known ones. Using the observed data on the fluxes of cosmic rays, it is then possible 
to constrain the DM annihilation rates. In this section, we employ the Fermi-LAT data on the 
mid-latitude (10° <| b\< 20°, 0° < / < 360°) y-rays IM and PAMELA data on the antiproton-to- 
proton flux ratio p/p f]\\\ to constrain our effective interactions. 
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A. Diffuse 7-rays and antiproton background estimation 



Gamma rays play an important role in indirect searches. As they are not deflected by galactic 
magnetic fields, they point back to their genuine sources. Several mechanisms may be responsible 
for the galactic diffuse gamma ray backgrounds: decays of mesons produced in nuclear in- 
teractions between cosmic rays and the nuclei in interstellar medium, inverse Compton scattering 
(IC) of electrons off photons and bremsstrahlung from electrons in the Coulomb field of nuclei 
[I72I1 . In addition to the galactic background, one also expects a contribution from the extragalactic 
background (EGRB). The sources of these gamma rays include other galaxies, unresolved point 
sources, large scale structures and interactions between ultra-high energy cosmic rays and CMB 
photons. Since each of them has different properties, it is difficult to predict the shape and magni- 
tude of EGRB. Observationally, the EGRB can be obtained by subtraction from the data. 

The positrons and antiprotons are also important indicaters of DM, because they are relatively 
rare in the galactic environment. As they traverse the interstellar space, they propagate randomly 
under the influence of galactic magnetic fields, and lose energy through processes of IC and syn- 
chrotron radiation [llOll . This process may be described by a complicated diffusion equation. The 
charged particles traversing the solar system are also affected by the solar wind, which results in 
a shift in the spectrum observed at the Earth compared to the interstellar one. For this, the solar 
modulation potential is taken as a free parameter ranging from 300 to 1000 MV 11731 17411. 

A good approach to estimate the spectra of the above cosmic rays is offered by the GALPROP 



code uM- It parameterizes the gas densities, nuclear cross sections and energy spectra for different 
processes, and then solves the diffusion equation numerically to get a complete solution for the 
density map of all primary and secondary nuclei. We have therefore used GALPROP to calculate 
both background and DM annihilation components in cosmic rays. 

For the galactic diffuse y-rays background, we calculate various contributions and fit the Fermi- 
LAT EGRB data with a power-law spectrum, E^d^/dE = 4>o(£'/GeV)^^. In calculating the 
antiproton galactic background, we scan the solar modulation potential to find out that the minimal 

is reached for the RAMELA p/p data at a value of 330 MV. The results are presented in Fig. U\ 
together with the corresponding to Fermi-LAT 7-rays and RAMELA p/p flux ratio data. 
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FIG. 7: Left panel: Background estimation for Fermi-LAT y-ray flux including various contributions and 
EGRB. The EGRB component is obtained by fitting Femii-LAT EGRB data with a power-low spectrum. 
The background flux (dark solid line) is obtained by adding all contributions. Right panel: Background 
estimation for R4MELA p/p flux ratio. The green dashed (orange short-dashed) line corresponds to the 
flux without solar modulation (with a solai^ modulation potential of 330 MV). 

B. Fermi-LAT and PAMELA bounds 



Now we derive indirect detection bounds using the Fermi-LAT mid-latitude y-rays and 
PAMELA p/p flux ratio data. As before, we treat one operator at a time. It turns out that each of 
the four operators contributes the same to the spin-summed and -averaged total cross section 
for the annihilation process — //, 




(16) 



Here s is the center-of-mass energy squared, rtif and M are respectively the masses of the final (/) 
and initial C¥) particles, and A^^ = 1 (3) when / is a lepton (quark). We have set one to unity 
and others to zero, is a dimensionless function. 
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(17) 



9 6 ■ 1087?2 1087?2 1087? 277? 54 ' 27 

with r = ruy/s and R = /s. From Eqs. (I16I17I) . one calculates the thermally averaged annihila 
tion cross section (c7|v|) to be 



1 ^ / mW 1 . o OS 50M'^-A9MWf+nm'^r ^ 



(18) 
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where v is the relative velocity of the annihilating DM particles, and the average is over the DM 
velocity distribution in the particular physical processes considered. 

The primary particles (quarks and leptons in our case) from DM annihilation will generate 
secondary particles (photons and antiprotons) at the production point via parton showers and 
hadronization. These processes can be simulated by Monte Carlo programs. For a given DM 
mass from 5 TeV to 10 TeV, we use PYTHIA6 .4 [76] to simulate dN^/dEy and dN^/dEp. Here 
dNy/dEy is the energy spectrum of photons produced per annihilation into the final state /, and 
dNp/dEp is that of antiprotons. We considered the — qq, channels for the photon spec- 
trum simulation, and — )■ qq for the antiproton. To ensure the accuracy, we took 10^ events in 
simulation. 



The photon flux in a given region AO. due to Dirac-type DM annihilation, is written as 
d^y _ 1 JAQ^, , dN; 



7ao dN'^ — f f 

4^iM2E(-M>/4. l=fj:i(>,A[j^pLM^,0)l (19) 

where / runs over all quark and lepton channels. r{s, 0) = (rg + 5^ — 2r05cos 0) is the galactic 
coordinate, r© the distance of the Sun to the galactic center, the angle between directions of 
observation and galactic center, and s the line of sight (l.o.s) distance. In terms of the galactic 
latitude b and longitude /, one has cos0 = cos cos/. We set the integral region in / to be the 
Fermi-LAT mid-latitude region (10° <| |< 20°, 0° < / < 360°). In our calculation, we assume 



the NFW profile [17 811: 

.n 2 

(20) 



Phalo('') rQ 



P© r 



l-fr/i? 



where pQ is the DM density at the solar location, and R the scale radius. We adopt the following 
values for these parameters: pQ = 0.3 GeV cm^^, = 8.33 kpc, and R = 20 kpc. The galactic 
DM particles should follow the Maxwell-Boltzmann velocity distribution. We choose the veloc- 
ity dispersion v = ^/ {v^{rQ)) = ^/3jlvc{rQ), with Vc{rQ) = 220 kms~^ being the local circular 
velocity, and thus (v^) = l^^/Jy^ifoj). For antiprotons, the source term is 

1 , ^ dNt 

= i^^'^halo('-)i;(c^|v|),^, (21) 

where q runs over all quark channels. We then implement Eqs. (1191211 ) into GALPROP to calculate 
the photon flux and p/p flux ratio with the same parameters in background estimation. 



To obtain the exclusion bounds, we adopt a simple statistical measurement in Ref. nJ% . Adding 



the DM component with the background flux to get the total flux, 4>totai = •I'bkgd + ^dm, we define 
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the difference t^^x^ = Ztotai(^' (^I^D) ~^min' where Ztotai(^' i^M)) the of the total flux and 
Xmin ^^^^ '•^^ background. The 90% CL limits are then obtained by requiring Ax^ = 2.71. Fig. [8] 
shows 90% CL limits on the thermally averaged cross section (c7|v|) from the Fermi-LAT y-rays 
and PAMELA p/p flux ratio data. For comparison, the upper limits from the ATLAS-SR4 and 
CMS-SR4 are also shown. 

There exist another y-ray limits based on the Fermi-LAT observation of Milky Way dwarf 



spheroidal satellite galaxies (dSphs) HSOll . By a joint likelihood analysis to 10 dSphs with 24 
months of Fermi-LAT data, those authors obtained the 95% CL upper limits on (Majorana-type) 
DM annihilation cross sections for several channels for a DM mass from 5 GeV to 1 TeV, assuming 
DM couples only to one channel at a time. They are about 10-^^-6.5 x lO^^s 

cm-^^s for the bb 

channel, 8.4 x 10^26 _ g ^ lO^^s^j^^Sg-i f^j. ^+^-^ and 1.4 x 10"^^ - lO^^^cm^s^^ for t+t". 
Upon multiplying them by a factor of two to convert to the case of Dirac-type DM, these limits 
should be compared to ours on the total cross section which has been obtained assuming that DM 
couples simultaneously and equally to all possible channels. Their most stringent limit in the bb 
channel is comparable to our PAMELA p/p limit, while their limits in the lepton channels 
and T+T^ are comparable to our Fermi-LAT limits but less stringent than our PAMELA p/p 



limit. One could also use the PAMELA positron fraction excess [18 ill to constrain the annihilation 



cross section. However, it has been suggested that the excess might arise from some astrophysical 



sources that were not accounted for earlier, leaving the origin of excess still unclear so far [|82 



VI. COMBINED CONSTRAINTS 

In the last sections we studied the individual constraints on DM interactions coming from LHC 
searches, direct and indirect detections. In this section, we put them together with the constraint 
from the DM relic density. 

The currently observed relic density is a remnant of DM production and annihilation in earlier 
epochs of our universe. If the annihilation was too fast, there would be not much DM left nowa- 
days; and in the opposite case, DM would be over dense in the current epoch. The observed value 
can therefore set a constraint on the DM interactions responsible for annihilation. We apply the 



standard procedure to calculate the relic density [|84i 18511 . i.e., by solving the Boltzmann equation 



numerically with the annihilation cross section in Eq. (|T6l ). summed over all fermions. 

In Fig. [9] we present the combined constraints for any one operator for DM mass from 
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FIG. 8 : 90% CL upper limits inferred by Fermi-LAT 7-rays and PAMELA pjp flux ratio on DM annihilation 
cross section as a function of mass M. For comparison, 90% CL limits by ATLAS and CMS experiments 
are also shown. 
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FIG. 9: Combined constraints on the effective scale A versus DM mass from 6 GeV to 1 TeV, including 
observed relic density, ATLAS experiment, SI XENONIOO direct detection, Fermi-LAT mid-latitude y- 
rays data and PAMELA pjp flux ratio data. The relic density bound is fixed by WMAP7-i-BAO-i-//o best-fit 
value, while all other bounds correspond to 90% lower limits. 
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6 GeV to 1 TeV. All limits correspond to 90% CL except that the relic density is fixed at 
WMAP7+BAO+//0 best-fit value, ^Duh^ = 0.1123 ±0.0035 im. With one species of DM, the 
latter would fix the effective scale A as a function of M. If instead one assumes that the particle 
under consideration is only one of the DM species in the universe, the relic density will set an 
upper bound on A. This is in contrast with all others which set lower limits. For light DM, LHC 
provides the most stringent constraint, while in the DM mass region from 30 GeV to 1 TeV, the SI 
direct detection by XENON 100 is most restrictive. 
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FIG. 10: Similai- to Fig.|9l but for DM mass from 6 GeV to 10 TeV. 



There is no direct detection constraint when the DM mass is larger than 1 TeV. However, one 
can still have collider and indirect detection constraints. In Fig.[lOl we show the combined con- 
straints from the relic density, LHC searches and indirect detections for DM mass from 6 GeV to 
10 TeV. The largest DM mass corresponding to the ATLAS limits is set to be 3 TeV. We find that 
in the region M > 1 TeV, the PAMELA pjp flux ratio is most stringent. In the LHC searches, 
as we discussed in secUlIl the cross section falls off rapidly in this mass region, resulting in very 
weak bounds on A. However, as the center of mass energy ^ increases, the fall-off will shift to 
even larger M. We therefore expect that LHC running at the projected energy of 14 TeV will be 
capable of pushing the bound on A further to higher values. 
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VII, CONCLUSION 



We have made a comprehensive analysis on the proposal that DM is composed of a spin-3/2 
particle which is a singlet of SM. Assuming certain kind of parity for such a particle or assuming 
the lepton and baryon numbers are still conserved, its leading effective interactions with ordinary 
particles would involve a pair of them and a pair of SM fermions. Demanding it to be a singlet turns 
out to be rather restrictive. There are only four types of effective operators in the form of products 
of chiral currents; and furthermore, all of them have the same or very close phenomenological 
effects thus simplifying significantly the physical analysis. 

Based on the above effective interactions we have investigated DM effects in various exper- 
iments and observations. These include the collider searches at LHC for monojet plus missing 
transverse energy events, direct detections by spin-independent and -dependent scattering off nu- 
clei, indirect detections via observations on y-rays and antiproton-to-proton flux ratio in cosmic 
rays, and the relic density. We found that the current data already set strong and complementary 
constraints. For a relatively light DM particle, say, below 30 GeV, where the DM pair production is 
much enhanced, the LHC experiments provide very stringent bounds. For instance, the latest data 
by ATLAS-7 TeV with an integrated luminosity of 4.7 fb~^ restricts the effective interaction scale 
to be above 15 TeV to 100 TeV for a dark matter mass of 20 GeV or so. The spin-independent 
detection by XENON 100, on the other hand, offers the most severe constraint for the mass range 
between 30 GeV and 1 TeV, where the effective scale is required to be above about 20 TeV. Rela- 
tively less severely constrained is the heavy DM scenario. For DM mass of 1-3 TeV, the strongest 
bound comes from the antiproton-to-proton flux ratio in cosmic rays, which excludes effective 
interactions with an effective scale lower than TeV. 
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Appendix: Amplitudes squared for subprocesses of dark matter pair production plus a monojet 



We list here the spin- and color-summed and -averaged amplitudes squared for the subprocesses 
studied in sec|llll They were used as the input file in our numerical simulation code. Denoting the 
momenta of the initial-state partons, the outgoing parton and the DM particles by pi^2, kj, and ki^2 
respectively, we define the kinematical variables 

2pa-kj = sxa, 2pfkt = sjb , (22) 

where 5 = (pi + P2Y, a, b assume values 1, 2, and j refers to jet. Ignoring parton masses and 
using kf2= M^, we can express other scalar products of momenta in terms of the variables, 

2pi-k2 = (1 -xi -yi)s, 2p2-ki = (1 -X2-y2)s, 

2ki-kj = {xi+yi-y2)s, 2k2-kj = (x2-yi+j2)*, 

2ki-k2 = {l-Xi-X2)s-2Nf. (23) 

For the sub-process q{pi)q{p2) — ^ ^(^/)^/i(^i)*?'v(^2), the amplitude squared is given in Eq. 
Q, where, for either of the operators ^1^2, B is given in Eq. (|7]), and for either of (^3^4, 



B 



m^{\ -xxf + 2M^s{\ -XI -X2)yi{2xi +lyi-2) 

-AM^s^{\-xi-x2fy'\ + s^{\-xi-x2fy]\+[xi ^ X2, yi ^y2). (24) 



For the sub-process q{pi)g{p2) '?(^7)^^(^i)^v(^2). we find 
where for ^1^2, 

B = 4M^[{l-xif + {xi+X2f]+2M^s{l-xi-X2){24xiyi-12xiy2 

-2x2yi + 2x2y2 + lOx? - 2^2X1 - IQxi + Uy\ + lyl - 12^1 - Uyiy2 + 5) 
-AM^s^{\ -XI -X2) 2 [Axiyi - 2xxy2 + 2x\ - 2xi+2y\+yl - 2yi - 2yxy2 + l) 
+5^ ( 1 - XI - ;C2) ^ [Axiyi - 2xiy2 + 2x\ - 2xi + 2^1 + - 2^1 - 2>'i3;2 + l) , (26) 
and for ^3,4, 

B = AM\{\-xif + {xi^X2f]-^2M\{\-xi-X2){2Axiyi-\2xiy2 

-2x2y\ + 2x2y2 + lO.^? - 2x2Xy - \<dxx + \Ay\^ly\ - \lyx - Uyxy2) 
-m^s^{\ -XX -X2f {4xiyi - 2xxy2 + 2x\ - 2xx + 2y\-^yl - 2yx - 2yxy2) 

W{\-xx- X2f {-2x2yx + 2x23^2 +^2 + ^yl- ^ym) ■ (27) 

20 



Finally, the result for the subprocess q{p2)g{pi) — )■ ^(fc;)^^(fci)^v(^2) can be obtained from Eq. 
(l25l) by the interchanges, xi ^ X2-, yi ^ yi- 
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